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ABSTRACT 



The state of the art in the development of the DPDS tech- 
nique at the Turbopropulsion Laboratory at the Naval Post- 
graduate School is reported. In principle, this method allows 
the determination of three dimensional flow fields behind 
rotors. A review of the instrumentation and hardware is given 
and two different data reduction procedures developed earlier 
are described. The results of reducing one set of rotor exit 
data using the two methods are compared. In order to verify 
the experimental procedures at low cost a bench-top periodic 
flow generator was used first to provide signals similar in 
type and frequency to those from the Transonic Compressor in 
which the rotor exit flow was measured. The data measured in 
the compressor were analyzed, and problems encountered in the 
routine application of the reduction methods were identified. 

The possible application of similar techniques in the Ro- 
tating Annular Cascade (RGG) at the DFVLR Gottingen is as- 
sessed and recommendations are made for the coordination of 
future experimental programs. 
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I. 



INTRODUCTION 



It is one important goal in today's axial turbine de- 
sign — for stationary steam and gas turbines in electricity 
generating plants as well as for jet engines — to develop ma- 
chines which are as compact as possible. Therefore larger 
enthalpy differences per stage are used in today's designs 
than was the case in the past. This results in supercritical 
pressure ratios for the conditions occurring in some parts of 
the blade rows, and leads to flows which are partly subsonic 
and partly supersonic. This regime is usually termed tran- 
sonic, and is well known to be difficult to investigate both 
theoretically and experimentally. This holds true for exter- 
nal flows — flow past single airfoils, for example, but it is 
particularly true for internal flows — flow within turbomachines 
being an extreme case. 

For large power-generating steam turbine plants the shaft 
speed is fixed by the RPM of the generator. Therefore an in- 
crease in work output is only possible by increasing the mass 
flow rate. As a result, progressively larger diameters of 
these kinds of machines are being used. Near the tip sections 
in the last stages of low pressure steam turbines this also 
results in transonic relative velocities. 

Presently, it is not possible to fully describe the pas- 
sage of the flow through the different stages of a turbine 
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from either experimental or theoretical approaches. In 
reality, very complicated three dimensional unsteady 
flow conditions exist. 

However, in order to simplify the picture for design pur- 
poses, two dimensional flow models are usually used. Experi- 
mentally, two dimensional cascades of axial blade row segments 
are generated either in an annular cascade of high hub-to-tip 
ratio or by developing the coaxial stream surface of interest 
into a plane surface as a rectilinear cascade. Because of 
the limited number of blades which must be used, the simula- 
tion of the required periodic flow field in a rectilinear 
cascade is usually imperfect, and this is particularly true 
for transonic flow conditions. 

This disadvantage is avoided in the case of an annular 
cascade, and for some years such a facility for Rotating 
Annular Cascades (RGG) has been in operation in the DFVLR- 
AVA in Gottingen, West Germany. In developing this facility, 
it was necessary to develop a new method of evaluating the 
measurements such that the results were representative of the 
flow through a truly two dimensional cascade. The method was 
required to be simple but reliable. Such a method was de- 
veloped and was verified by investigations on two cascades of 
very different types [1, 2]. 

The fundamental advantage of the method is that the ex- 
perimental determination of only the total pressure and 
total temperature in the homogenous downstream flow field is 
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necessary. This requires the existence in the downstream 
flow field of a region where non-uniformities from wakes have 
disappeared and where the boundary layers from hub and tip of 
the annulus have not penetrated. Although the existence of 
such a region was shown it was considered desirable to mea- 
sure peripheral and time-varying flow properties as well. 

In order to measure the unsteady flow field behind a 
rotor one must either use high response probes or an optical 
technique of some kind. While laser velocimetry techniques 
have many advantages, they are normally able to measure only 
two components of the velocity. They do not provide measure- 
ments of the pressure field, and in highly unsteady flows the 
accuracy must be questioned. 

At the Turbopropulsion Laboratory (TPL) of the Naval 
Postgraduate School (NPS) , the dual-probe digital sampling 
(DPDS) technique is being developed. The technique is based 
on the use of two Kulite semi-conductor pressure probes of 
very simple design. The approach has the potential of evalu- 
ating all three components of the periodic velocity vector in 
a turbomachine or rotating cascade. Following early work 
with a single probe in a compressor flow field [21], effort 
at NPS was devoted to developing and verifying a method of 
calibrating and applying the dual-probe system in the steady 
flow of a free jet [3]. The application of the probe system 
in the compressor flow field was delayed until recently by 
drive turbine bearing failure. 
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Now, some first data from the compressor and early at- 
tempts to analyze the data are reported here. While there are 
some open questions left by the work done so far, it is recom- 
mended nevertheless that the method be applied in the RGG at 
the DFVLR. Because of the early stage of development and 
complexity involved in the process it is suggested that the 
efforts at the DFVLR and NPS be conducted in close association 
so that improvements are directly available to both institutions. 

This account should be viewed as a documentation of the 
learning process that the author has undergone. It reviews 
the earlier work listed as references. Key passages from 
those references have been included without change on occa- 
sion when any contraction of the wording would have detracted 
from the descriptions originally given. The list of references 
is organized so as to group the materials logically, rather 
than in the order that they are mentioned in the text. The 
references are reviewed in the following section. In sections 
III and IV, descriptions of both the experimental and computer 
hardware are given. Section V describes the DPDS-method and 
in section VI some results are shown. Finally, in section VII 
ideas for a first application in the RGG are described. 
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II. REVIEW OF PUBLICATIONS DEALING WITH 



UNSTEADY FLOW MEASUREMENTS IN TURBOMACHINES 

In [3] to [43] publications are listed which are a useful 
source of information in the field of unsteady flow measure- 
ments in turbomachines. The quoted references are those that 
the author has used; the list does not cover the available 
literature completely. It should be mentioned that in the 
following text, the references referred to most often are 
those which concern the research activities of the TPL. 

The sequence of the reference list is such that Masters 
Theses of the NPS ([4] to [10]) are given first, followed by 
Technical Reports of the NPS ([12] to [17]) and finally, 
publications in the open literature ([18] to [21] and [3]), 
which generally summarize much of the useful material from 
the theses and reports. In [11] a U.S. patent is referred 
to. References from [22] to [43] deal mostly with different 
kinds of high response probes used in the flow field down- 
stream of rotors. Some results from laser measurements are 
reported in [37] and the use of Schlieren pictures for in- 
terpreting the flow field in [40]. In [31], [39], [41] and 

[42] the use of high response pressure transducers to measure 
the dynamic flow behavior on rotor blades and stator blades 
are reported. Case wall measurements based on the same kind 
of instrumentation are reported in [36] and in [32] a FORTRAN 
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program is given for reducing case wall static pressure data 
to contour maps. In [24], a general description of different 
problems which can occur in investigating flow fields in 
transonic compressors, is given. 

The ability to measure the true time-averaged static and 
impact pressures is of considerable importance and is needed 
to permit a full interpretation of semi-conductor probe 
measurements in non-steady flow fields. This is discussed 
in [26], [27], [33] and [34]. For application in transonic 

flow fields the oil-filled total pressure probe described 
in [33] and [34] may be the most useful tool available. 
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Ill . REVIEW OF THE EXPERIMENTAL EQUIPMENT 



A. TRANSONIC COMPRESSOR (TX) 

The transonic compressor is located in the second test 
cell of the High Speed Turbine and Compressor Laboratory. 

This laboratory in total consists of three test cells, a 
control and instrumentation room and a compressor room con- 
taining a single 12-stage axial air compressor that provides 
the power for all experiments. ’ There is also a "spin pit" 
for the determination of stresses in rotors. Details of the 
laboratory are reported in [44]. 

Compressed air at up to 11 lbs. per sec. at three atmos- 
pheres pressure can be valved to each cell from an Allis- 
Chalmers multi-stage axial compressor. The 1250 HP compressor 
has a variable speed drive with an automatic control and a 
by-pass surge suppression device. The mass flow rate from 
the compressor is metered by a sharp edged orifice meter lo- 
cated downstream of an aftercooler in the outlet pipe. 

A line drawing of the TX is shown in Fig. 1 . The compres- 
sor is a single stage axial fan designed to operate at 30,460 
RPM with a relative tip Mach number of 1.5. The rotor diameter 
is 27.94 cms (11 inches) , and the hub-to-tip ratio at the 
rotor face is 0.5. At the design RPM and referred flow rate 
of 8.62 kg/sec (19 Ib/sec) , the relative flow angle at the 
rotor tip is 65° and the design pressure ratio is 1.6. The 
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compressor rig is open loop with an air filter, throttle and 
flow nozzle in the induction pipe and with the compressor 
exhaust going to atmosphere via a radial diffuser. The flow 
rate is controlled by an electro-hydraulic rotating throttle 
plate located in the inlet duct. The inlet of the transonic 
compressor is designed to contract and to rapidly accelerate 
the flow in order to generate near uniform conditions, with 
thin boundary layers, at the rotor face. 

The compressor test rig has the advantage of providing 
redundancy in measurements. The operating power can be de- 
termined in three ways: either from the mass flow rate and 

temperature drop through the turbine, or from the mass flow 
rate and temperature increase through the compressor, or from 
the speed and the torque measured by strain gauges on flex- 
ures between the compressor stator (mounted on bearings) and 
the machine casing. 

B. INSTRUMENTATION 

1 . Probes for Radial Surveys 

Probe survey ports are provided at four axial planes 
and fixed Kiel and temperature probes are installed for over- 
all stage performance measurements as shown in Fig. 1. The 
case of the compressor can be rotated manually about the axis 
while the machine is running in order to position probes and 
sensors mounted in the case wall so as to obtain surveys in 
the peripheral direction (with respect to the fixed stator 
blading) . 



17 



The combination temperature-pneumatic probe which is 



used is shown in Fig. 2 . One of those probes is used in front 
of the rotor and the other one is located between the rotor 
and the stator. In use, the probes must be rotated about 
their axis to balance the pressures at the side sensors, P 2 
and p^, and the yaw angle can be read from a vernier scale. 

The center sensor, p^^, measures an indicated impact pressure, 
the pressures P 2 and p^ approach static pressure and the 
fourth sensor, p^, provides a measurement of the pitch angle, 
when the probe has been suitably calibrated. 

The probe was designed specifically to be used down- 
stream of the transonic rotor [4]. The use of separate sen- 
sors (the "claw" arrangement rather than pressure taps on a 
shaped surface) was to minimize the "immersion" effect. It 
is also likely that the separate sensors would be less in- 
fluenced by unsteady effects from the rotor. The separation 
of the sensors peripherally is of no consequence because sen- 
sors at the same radius would automatically average in the 
same passing flow field. The pitch angles in the compressor 
are expected to be small and severe velocity gradients in the 
radial direction are not expected away from the hub and tip 
regions. The necessary radial separation of the fourth sen- 
sor is therefore not considered to be critical. Four sensors 
are used rather than five, to allow the simple arrangement of 
tubes shown in Fig. 2. This arrangement minimized blockage 
in a strong, simple, leak-free construction and gave suitably 
short response times. 
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The pneumatic pressures are sensed by a conventional 
differential strain gauge transducer and Scanivalve referenced 
to atmosphere. 

The separation of the thermocouple sensor from the im- 
pact pressure sensor in the radial direction. Fig. 2, is also 
not important since it could be properly accounted for in the 
data reduction from radial surveys. 

Probe calibration procedure and results, data from ra- 
dial surveys and stage performance measurements are reported 
in [17]. 

In [21] a method was described for using two semi- 
conductor pressure probes together with the technique of syn- 
chronized sampling, to obtain the distribution of the velocity 
vector downstream of a rotor. The geometries of the two probes, 
designated Type A and Type B, and their installation in the 
compressor annulus are shown in Fig. 3 . The two probes are 
positioned as shown at the same radial displacement, but sepa- 
rated peripherally around the machine annulus. The probes are 
mounted in such a way that they can be rotated about an axis 
through the probe tip. The details of the geometry of the two 
probes are shown in Fig. 4 . 

The particular shape of the 3 probe is chosen to enable 
measurements to be made close to the case wall, and therefore 
limitations occur in the measurements which can be made close 
to the hub . 
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It was argued that, in principle, by rotating the 
probes in yaw about their tips and controlling the sampling 
of the data from each probe to be at the same position in the 
rotor frame, the system of two separate probes could be used 
to acquire data at a single point in the periodic flow from 
the rotor, thus obtaining data corresponding to those obtained 
from the sensors of four- or five-hole pneumatic probes when 
measuring velocity in steady flows. The technique involved 
the use of probes having the simplest geometry and thus avoided 
the large size, expense and unreliability of multiple sensor 
probes which incorporate multiple semiconductor transducers. 

It would be far cheaper to replace a single probe of simple 
geometry in the event of a failure compared to repairing or 
replacing a multiple sensor probe. Because of the simple sen- 
sor tip geometry (that of a cylinder at incidence to the flow) 
the unsteady response was likely to be as good as could be ex- 
pected of any single physical sensor. 

All probes are yawed and moved to different radial 
positions, manually. Both yaw angle and radial position are' 
recorded automatically using position potentiometer circuits. 

2 . Wall Pressure Measurement 

Twelve Kulite CQL-080-25 pressure transducers with 
natural frequency about 125 kHz are mounted with their dia- 
phragms flush with the inner case wall of the compressor. 

The transducers are located along four axial lines separated 
peripherally by 10 and 90 degree arcs as shown in Fig. 5. 
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The transducers are used in conjunction with Datel Model 20 1C 
instrumentation amplifiers which have a flat frequency response 
to 100 kHz. Each Kulite pressure transducer is matched by a 
pneumatic static pressure tap at the same axial location in 
the case wall, but displaced circiomf erentially . Other pneu- 
matic static pressure taps are provided upstream and downstream. 
The reference tubes of the transducers are connected by separate 
plastic tubes to a manifold fed from a controlled pressure 
source or to the pneumatic line from the pressure tap at the 
same axial location. Clamps were used to valve the reference 
side of the transducer either to the regulated pressure or to 
the corresponding pneumatic pressure tap. 

3 . The Timing Disk 

To enable synchronization of the sampling of the semi- 
conductor pressure transducer outputs with the rotation of the 
rotor, a timing disk is fitted to the rotor shaft. Fig. 6 . The 
disk has a circle of eighteen equally spaced holes (equal to 
the number of rotor blades) at one radius and a single hole at 
a smaller radius. Light emitting diodes and photo detectors 
aligned with the holes generate bell-shaped electrical pulses 
at one-per-revolution and one-per-blade frequencies. The sys- 
tem is described in detail in [5]. 

The blade passing frequency, f, can be computed by 
multiplying the wheel speed (revolutions per second) with the 
number of blades. The maximum rpm of the TX is about 500/sec 
and the rotor constsis of 18 blades. Therefore, the e.xpected 
frequency range is up to 10 kHz. 
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C. DEVICES FOR DATA ACQUISITION 



The primary hardware devices used for data acquisition are 
the HP-21MX computer with its 20 Megabyte magnetic disc, print- 
er, plotter, and two terminals, the "Pacer" and the HP-5610A 
16-channel A/D converter. The complete data acquisition sys- 
tem is shown in Fig. 7. 

1 . Hewlett Packard HP-21MX Computer 

The HP-21MX is a (Micro-programmable) mini-computer 
having 128 machine instructions and 32K of logical in 24QK of 
main frame memory. A 20 megabyte capacity disc and disc drive 
surmount the core memory limitation and provide multiple use 
under time sharing. 

An important feature which is typical of computers of 
this size is the input-output structure. With a limited num- 
ber of relatively slow I/O devices to be serviced, the com- 
puter can communicate with all devices through a single port 
known as the I/O bus. Each device requires its own I/O inter- 
face on the bus. The interface acts as a filter and ensures 
that output information is received only by the device desig- 
nated to receive it and that input information is put on the 
bus from only one device at a time. The I/O software drivers 
control the I/O hardware interfaces by commands to either 
"turn-on" or "turn-off". 

2 . Pacer 

This peripheral device which v;as designed and con- 
structed to trigger data acquisition from a stationary sensor 
at any fixed point in the rotating rotor frame, independent 
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of the rotor speed, is shown schematically in Fig. 8 . In its 
original form, a detailed description of the internal opera- 
tion is given in [6]. 

The bell-shaped optical signals from the timing disk 
are first converted in the "wave shaper" section into square 
pulses suitable for TTL (Transisi tor-Trans is tor-Logic) circuits. 

The "Pacer" consists of two major sections, a "synchro- 
nized command pulse section" and a "RPM counting section". 

The "synchronized command pulse" section is the heart 
of the "Pacer" . The one-per-rev. signal is used as a master 
pulse in the timing process. The one-per-blade signal is 
halved in frequency through a counter and input into a PLL 
(Phase-Locked Loop) circuit. The PLL "locks-on" and outputs 
a similar wave form at 256 times the input frequency. A com- 
parator compares the number of counts following the master re- 
set (taken from the one-per-rev. signal) with the programmed 
number received from the computer. When they agree, the com- 
mand signal from the computer to the A/D converter is allowed 
to pass and conversion occurs. 

The device can be programmed to time the data conver- 
sion to occur at any point in the rotor frame, since the "lock- 
on" makes it independent of the rotor speed. 

The pacer can also be programmed to allow conversion 
either at a single location in a particular blade passage, or 
at the corresponding positions in each blade passage. 

The "RPM counting section" continuously counts the num- 
ber of 250 kHz time base pulses (obtained from the computer's 
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1 MHz clock) that occur between the one-per-revolution pulses 
received from the test rig. This number of counts is available 
as an output from the "Pacer" on every revolution of the rotor. 

The original "Pacer" was recently modified to overcome 
the following three limitations which were realized when using 
the original version: 

(1) The timing pulses generated with the "Pacer" were 
not always spaced linearly in time between blade pair synchro- 
nizing pulses. This resulted in data which in some cases was 
subtly distorted, and in other cases appeared to have noise 
riding on it. 

(2) The range over which the "Pacer" could follow 
rotor RPM changes and remain synchronized was limited to ap- 
proximately ±15% of the initial RPM at which the "Pacer" was 
set to take the data. This required a careful manual adjust- 
ment of trimming capacitors and resistors to obtain "lock-on" 
in each range of RPM in order to acquire data at different 
speeds . 

(3) The rate at which data could be taken was limited. 
Rather than being able to sample data on every revolution of 
the rotor, the system was only capable of taking data once 
every 3 to 10 revolutions, depending on RPM. 

Changes to the "Pacer" are described in detail in [10] 
which fully overcame the above three limitations. They in- 
volved changes in both hardware and software over those re- 
ported in [6]. The present system is extremely easy to use 
and presents no obvious limitations. 



24 



3 . Hewlett Packard HP-5610A A/D Converter 



Up to sixteen channels of analog input voltages in the 
range of ±1.0 volt can be sampled and recorded by the A/D con- 
verter under the control of the 21MX computer. The acquisition 
process is governed by the built-in features of the A/D con- 
verter and by the programs which can be written for the 21MX. 

The HP-5610A Analog-to-Digital Converter is designed 
to convert analog voltages into ten data BIT's at the maximum 
throughput rate of 100 kHz. The HP-5610A can operate in one 
of six modes. The two fundamental modes of operation are "se- 
quential" — in which all 16 channels are scanned in numerical 
order, and "random access", in which a specific channel is 
sampled on receipt of a command word and an encode command 
pulse from the 21MX computer. 

Variations within these two basic modes are also pos- 
sible which allow for externally supplied sequencing pulses 
and for externally controlled timing of samples. 

Another mode which is used only for non-paced data is 
the "Free Run, Random .Access" mode. In this mode the command 
word is required as before, but no encode command is needed 
from the computer. The A/D converter simply converts data as 
fast as it can (100,000 samples per second) on the selected 
channel . 

The command word tells the A/D converter which mode of 
operation to use and which channel number to sample. The en- 
code command pulse triggers the data conversion to start 2 tisec 
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later. The data conversion itself is finished in a total time 



of IQ ysec. Using computer-issued encodes, which is the mode 
required for paced data, the sample cycle time is 20 ysec. 
Hence data can be converted at rates of up to 5Q,QQ0 samples 
per second, depending on how rapidly each successive command 
word is received. With an input conversion aperture of 5Q 
nanoseconds, rapidly changing signals CIQQ kHzl can be con- 
verted accurately. 

The A/D digitizes data using a method of successive 
approximations. That is, a known test voltage is compared 
with the sample voltage beginning with 1/2 full scale (Q.5v 
in the present system) and a decision is made whether the test 
voltage is high or low. The comparisons continue, halving the 
range each time until the machine accuracy is reached. For a 
resolution of ±1.0 volt to IQ BIT ' s the accuracy of the A/D 
converter is 1/512 = 1.95 mV. 
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IV. DUAL-PROBE DIGITAL SAMPLING (DPDS) TECHNIQUE 



Until the development of miniature semiconductor trans- 
ducers having large bandwidth frequency response, real-time 
pressure measurements in high speed compressors were not 
practicable. Currently, using these pressure transducers, 
coupled with fast analog-to-digital conversion and computer 
controlled time delay circuits, a variety of new real-time 
measurement techniques are possible. The well-known tempera- 
ture dependence of the transducers (see for example [39]), 
can be overcome by an on-line calibration (OLC) such as is 
explained in Section IV C. 

The difficulties to be overcome in the design of a sam- 
pling method include achieving timing accuracy at high blade 
passing frequencies and accounting for the small variations 
in rotor speed that are always present. 

The technique of DPDS applied in the transonic compressor 
is illustrated in Fig. 6 and Fig. 9 . Referring to these 
figures, a transducer which is positioned downstream of the 
rotor and fixed to the machine frame outputs a time varying 
signal which is basically periodic at the blade-passing fre- 
quency. By successively triggering the sampling of the trans- 
ducer with a fixed delay following a one-per-rev. trigger 
pulse, the condition at a fixed point in the rotor frame can 
be determined. By programming the delay time, in principle 
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any chosen location in the rotor frame can be examined. How- 
ever, at design speed, the blade passing frequency is near 10 
kHz, and conventional time delay techniques, in the presence 
of speed fluctuations, can not provide the required resolution. 
A computer peripheral device called "Pacer" was therefore de- 
signed and developed to provide the programmable delay ([6] 
and [ 10 ] ) . 

If the stagnation temperature is known, four quantities 
are sufficient to determine the velocity vector: for example, 

static pressure, p^, total pressure, p^, yaw angle, a, and 
pitch angle, Following calibration however, any four inde- 

pendent measurements of pressures can be used to evaluate these 
same unknowns (as is done typically with a multiple sensor 
probe in steady flow [17]). In principle, the required mea- 
surements can be obtained by rotating the probes about their 
tips to four different directions where pressure readings can 
be taken. In reality, one must pay attention to the character- 
istic behavior of the probes, i.e. the output of the probes as 
a function of yaw angle, pitch angle and Mach number (M) in 
the ranges of these parameters which are to be measured in the 
compressor flow field. Experience shows in fact that it is 
preferable to evaluate curves of transducer pressure (p) as a 

function of probe yaw angle setting (a ) and to refer these to 

P 

calibration curves. It is obvious that probe readings at only 
four different yaw positions can not provide a precisely de- 
fined curve and therefore a comparison which is sufficiently 



28 



accurate. Certainly, more than four different probe settings 
are necessary. Seven positions were suggested in [3]; nine 
and eleven, respectively, are used for data reported here. 

In principle, there are many possible methods for reducing 
data acquired in the DPDS technique. To date three entirely 
different approaches have been attempted at the TPL. The 
first. Method I, uses a particular set of polynomial approxi- 
mations for the characteristics of the two probes obtained by 
calibration [3]. The second. Method II, is described in [14], 
and requires the representation of the calibration as an array 
of pressure coefficients. The method is a further development 
of the work reported in [13] and overcomes some of the earlier 
limitations. The description and results of the third method. 
Method III, which makes use of work reported in [15] in an ap- 
proach which parallels that used for steady-state velocity 
probes [17], is given in [45]. 

In this report, TX data are analyzed using Method I and 
Method II only. Some details of the procedures will be given. 



A. METHOD I — ANALYTIC.AL DATA REDUCTION [3] 



Using the probe pressure coefficient. 



(c ) , defined as 
P 



c 

P 




( 1 ) 



where p = pressure at the probe sensor, p^ = static pressure 
and p^ = total pressure, the characteristic behavior of the 
A-probe as a result of observation (first) and calibration 



29 



(second) can be analytically expressed by the following three 
fourth order polynomial equations; 
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where is the flow yaw angle with respect to the axis of 
the probe, a ^ indicates the value of where static pres- 
sure is measured and (p is the pitch angle. *^pq^ defines the 

maximum c value. 

P 

For the B-probe, only the equivalent to equation (3) is 
useful because its inclination (see Fig. 4) is such that the 
angle at which static pressure is registered, does not 

exist at most positive angles of pitch. The fourth expression 
is therefore 

CpoB “ ^0 ^l”^ (5) 

Af, C^, and (i = 0 to 4) are polynomial coefficients 
which are established by calibration. In Section V some more 
details are given when the evaluation of these curves for two 
new A- and B-probes is described. 

In the above equations for every probe yaw angle position 
which is set in the data acquisition procedure (up to 11 
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positions have been used) a single value for each of a, 4), 
and p^ can be derived. 

The flow yaw angle, a, is taken simply as the position of 
the maximum of the probe output when the probe pressure is 
plotted versus the probe yaw angle setting in the machine, 

(m) . Of course, this could be done for the output of both 
probes and, in the ideal case, the same value would be cal- 
culated. To apply the calibration equations the yaw angles 
must first be calculated relative to the axis of the probe; 
i.e. (m) = ap^(m) - a. 

The pitch angle, 4> » can then be calculated iteratively for 
every data point, m, from the following implicit equation for 
4> obtained using Eq. (1) and Eq. (2): 
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where Cp^^ and CpQ^ are given by Eq. (3) and Eq . (5) respec- 
tively, and is given by Eq. (4). p^^ and p^^ 

maximum values evaluated by fourth order curve fits through 
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the m data points, for the type A and type B probes respec- 
tively. It should be noted that only the maximum value of the 
B-probe is used in this procedure. The solution of Eq. (6) to 
obtain (}>, can be carried out using Newtonian iteration. 

The static pressure for every data point, p (m) , can be 

o 

calculated using 



, , POA - Pob'*'*’ 

p (m) = 

1 - 



( 8 ) 



By applying Eq. (1) to Pq^ the total pressure for every data 
point, p^(m), can be obtained using 



p^(m) = + 



Pqa - Ps^”* 

'^POA 



(9) 



It can be seen from Eq. (6) that for p^(m) = Pq^ the left 
side of the equation becomes unity. In the ideal case (per- 
fectly symmetric probe tip geometry) a^(m) should be equal to 
0 when p^(m) = Aq = 1, and no solution for ({) (m) is then 

possible. 

B. METHOD II — NUMERICAL DATA REDUCTION [14] 

Method II uses the values of c as defined by Eq. (1) and 

P 

obtained from calibration data in the form of arrays corre- 
sponding to a particular grid of yaw and pitch angles. Since 
stem effects and manufacturing imperfections will lead to de- 
partures from the perfectly symmetric probe behavior in yaw, 
a method which does not require an analytical expression for 
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the probe characteristics has potential advantages. The 
authors of [14] developed a numerical solution which avoided 
any requirement for curve fitting and, in principle, allowed 
a determination of p^, p^, a and <}) with only four measure- 
ments. In the current version, the pressure coefficients 
defined in Eq. (1) are assumed to be totally insensitive to 
Mach number and Reynolds number. Nevertheless, if any sig- 
nificant variations did occur an additional iteration scheme 
could be added to the computer solution. 

From four different measurements — three from the A-probe 
and one from the B-probe are suggested — four c^ values can be 
calculated. It is assumed that each characteristic has only 
one maximum in the range of interest and therefore one com- 
bination of yaw and pitch angle Ca, ;])) exists where the four 

c values fit those in the calibration tables. The FORTRAN 
P 

computer program given in [14] provides a systematic trial- 

and-error search process, essentially a convergence scheme 

on the two variables yaw angle and pitch angle. A linear, 

double-interpolation scheme is available as well as a second- 

order, double-interpolation scheme. 

The computer program in [14] was entered into the HP-21MX 

computer together with the c data arrays, also given in [14], 

P 

which are based on analytical functions. The c^ values so 
derived are automatically symmetric for positive and negative 
yaw angles, a. 



First, the program was run for the examples given in [14] 
The results agreed perfectly. Second, the method was applied 
using the actual (unsymmetric) c calibration data arrays. 

Ir 

Incomprehensible and obviously incorrect results occurred. 

After changing row and column indices of the c arrays in the 

P 

FORTRAN program (referring to [14], statements 3010 to 3040, 
page 31 and statements 4830 to 4900, page 44; in addition, 
IFLAG = 0 is omitted after statement 4990) results were 
evaluated, which were in the expected range. On rerunning 
the corrected program using the example data, only minor de- 
viations from the first results were found. Hence the cor- 
rected program is thought to be working for non-symmetric 
probe characteristics. Additional checks using calibration 
data from two new probes are reported in Section V. 

C. ON-LINE CALIBRATION (OLC) 

The reduction procedures described in the previous two 
sections require the input of correct values of the probe 
pressures. Because of the temperature sensitivity of the 
transducers, which affects mainly the zero balance but to a 
lesser degree the output slope of the strain gauge, correc- 
tions have to be made to the transducer voltages to obtain 
reliable quantitative pressure measurements. The problem 
occurs even though a temperature compensation module is in- 
cluded in the transducer circuit. The problem of temperature 
dependence is widely discussed in the literature; for example 
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[8], [23], [33], [39], Two approaches are in common use; 

either the signals are used as a measure of the fluctuations 
only and the time-averaged pressure level is measured by a 
second device, or flow temperature is measured simultaneously 
and corrections are computed based on an earlier temperature 
calibration. 

Progressive improvements have been made by the transducer 
manufacturers to reduce temperature sensitivity during the 
past decade. Also, in one case, automatic electronic com- 
pensation similar to linearization of hot wire circuits, is 
reported [ 46 ] . 

In the TX an on-line calibration procedure for all Kulite 
transducers (A- and B-probe and wall statics) is used to ac- 
count for temperature dependent effects. For the A-probe the 
method is described in detail in [8] and, therefore, only a 
brief summary is given here. 

To evaluate the slope (differential pressure vs. voltage 
output) , x^, of the transducers in situ, while ready to take 
data, the transducers are referenced to different pressures 
(6 at maximum) and a linear least squares routine is used to 
obtain x^. Experience shows that almost without exception 
the slope is linear and reproducible to within the accuracy 
of measurement for the modest changes in temperature involved 
in the experiments (''10-20°0 . 

Additionally, the above procedure establishes an inter- 
cept value, Xq, which can be monitored over longer times to 
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examine any temperature drift which may occur while test con- 
ditions are held constant. Experience shows that normally 
deviations equivalent to less than .5 in H 2 O pressure occur 
while dynamic pressures are never less than 20 in H 2 O. Larger 
deviations can be interpreted as poor temperature compensa- 
tion, which was the case, on one occasion, for one of the case 
wall transducers which was subsequently replaced. 

In order to interpret the intercept, Xq, the assumption 
is made that a pneumatic probe measures the correct time 
averaged pressure when positioned at the same radial location 
and set to the same yaw angle as the high response probe. 

This is to be regarded as a delicate problem because investi- 
gations are published, for example [33], [34], which clearly 

indicate that this can be a totally wrong assumption. To 
date, it has been of minor importance for the TX measurements 
due to low peak-to-peak pressure amplitudes which have been 
measured so far. An example is given in [17] and there it is 
indicated that the parameter "ratio of wave amplitude to mean 
value," on which the error depends, is 10 times less than the 
lowest value investigated in [34]. For future measurements 
at higher rotational speeds, care will be taken to reexamine 
the importance of averaging errors under the then prevailing 
machine conditions. 

When the transducers are referenced to the ti.me averaged 
pneumatic pressure their outputs are set to "0". In applying 
the calibration, samples of Kulite output voltage are reduced 
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to values of pressure, Ppj-Qj^g^ using the calibration coeff- 
cients and x^ in the equation 

Pprobe = ■'o ® Pref 

where Xq is expected to be very close to zero. In acquiring 
data, the transducer reference tube is usually open to atmos- 
phere and the reference pressure, therefore equals at- 

mospheric pressure. 

The total pressure, p^, read by the combination probe. 

Fig. 2, is used for the A-probe to evaluate the intercept, 

Xq . The intercept for the B-probe can not be found in this 
way because of its inclination. Fig. 4. Rather, the static 
pressure, p^, and flow pitch angle, b , from the combination 
probe are used in Eq. (5) to compute the maximum output of 
the B-probe, Pqq* This procedure is necessary to adjust the 
transducer of the B-probe to "0" on-line. It is noted that 
the process does not involve equating to a pressure which is 
measured directly in the TX. In principle a pneumatic equi- 
valent to the B-probe geometry could be used, but it is feared 
that minor differences of the inclination angle or tip geome- 
try could lead to totally unacceptable errors. The present 
procedure is one which ensures that the time-averaged results 
of the Type A and B probe system are the same as those given 
by the combination probe. Within this constraint, the Type A 
and 3 system can resolve the time-dependent behavior as a 
blade- to-blade distribution of flow properties. 



V. SOME RESULTS FROM HIGH RESPONSE PRESSURE PROBES 



A. PRELIMINARY INVESTIGATIONS 

1. Results from Periodic Flow Generator (PFG) 

A small periodic flow generator (PFG) was used by 
the author as a means of becoming acquainted with the high 
speed data acquisition system in the simplest way. The ap- 
paratus is described in detail in [5] and [7] and only the 
main features are given here. 

Regulated air is delivered through a settling cham- 
ber to a small nozzle. A sprocket (or interrupter) wheel 
driven by a speed-regulated electric motor, is mounted di- 
rectly above the jet such that the teeth of the sprocket pass 
through the jet. The spaces between the teeth simulate blade 
passages and "one per revolution" and "one per blade" signals 
are generated from a similar arrangement of holes, light 
emitters and detectors as is used on the TX. With 60 teeth 
on the wheel and 10,000 rpra from the motor, blade passing 
frequencies up to 10 kHz are obtained to simulate the condi- 
tions expected in the TX. The high response pressure trans- 
ducer is mounted on the centerline of the jet behind the 
plane of the sprocket wheels, so that, ignoring transient 
flow deflections in the unsteady process, the jet impact 
pressure is measured. In the potential core of the free 
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jet, the transducer would be expected to register a pressure 
equal to the supply pressure. 

In addition to tests conducted solely for the pur- 
pose of learning, several small experiments were carried 
out of some general interest. For example, results were 
obtained to document how well the "Pacer" works. In Figs . 

10 and 11 two examples are shown which are results obtained 
to verify that the "Pacer" was performing correctly as a 
trigger device. On the interrupter wheel one "blade" pair 
was distinguished by fixing a wire across the gap between 
two teeth. 3y using the "Pacer" to obtain measurements be- 
hind the different "blade" pairs, the blade space contain- 
ing the wire would be detected by the appearance of a 
"wake". From the two figures, which show the results for 
two different frequencies, it can be seen that the wake was 
present only for pair 2. All the other signals were quite 
similar, as also were those for pairs 5 to 9 which are not 
shown. 

At this stage it was not possible to take samples 
every revolution and, therefore, the test time was much 
longer than expected. This problem has since been solved 
[10] and examples are given below when measurements from 
the TX are discussed. 

In Fig. 12 results of taking paced (synchronized) 
samples of a sine-wave, produced by a signal generator, are 
given for four different numbers of repetitions. Again, the 
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"Pacer" was used as the trigger device. The number of re- 
petitions simply indicates how often a sample at a particu- 
lar location was taken. An average value was calculated 
and is shown plotted in Fig. 12. A trace of non-uniformity 
is seen only for the case of a single sample. 

The PFG was also used to obtain comparisons of the 
response of different kinds of pnevimatic and high response 
probes similar to those given in [7]. The goal was to get 
additional information to add to the work reported in [26], 
[33] and [34], for example, particularly to obtain results 
at higher frequencies. Unfortunately, the small dimensions 
of the jet and wheel did not allow controlled variations of 
frequency alone to be made and, again, the learning experi- 
ence was the main benefit derived from the experiment. In 
Section VI suggestions are given on how more detailed in- 
vestigations could possibly be carried out in the RGG at 
the DFVLR. 

2 . Calibration of New Type A and Type B Probes 

In [3], the procedure for calibrating the Type A 
and Type B probes in a free jet is described. The old probes 
in [3] were the ones shown in Fig. 4. Two new probes were 
bought which were smaller in size in the sensing area. The 
diameter of the head was 1.63 mm Cinstead of 2.36 mm). Addi- 
tionally, the inclination angle of the B probe was 35 degrees 
(instead of 55 degrees) as was used in [3]. The new probes 
were both derived from Xulite Type XB-062-25 transducers. 
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The probes were calibrated in the 107 mm diameter 



free jet exhausting into the room and supplied from the 
laboratory compressed air system. A sketch of the calibra- 
tion apparatus is shown in Fig. 13 . The probes were mounted 
in a probe holder which allowed adjustments to be made in 
yaw and pitch of the probes about the probe tips. The Mach 
number of the jet was controlled from approximately .2 to .7 
in steps of .1. The Mach number was obtained from the stag- 
nation pressure measured just upstream of the freejet nozzle, 
and atmospheric static pressure. The upstream stagnation 
pressure was determined to be equal to the stagnation pres- 
sure at the test probe station to within measurement accura- 
cy. All pneumatic pressures were recorded using a single 
strain gauge transducer, signal conditioning and Scanivalve, 
with the Steady State Data Acquisition System shown in the 
upper half of Fig. 7. The reference sides of the Kulite 
transducers were vented to atmosphere or supplied with con- 
trolled air pressure (also recorded through the Scanivalve 
system) . The yaw angle of the probes was recorded by the 
data system from rotational-position potentiometers to an 
uncertainty of about 0 -. 0 5°. The pitch angles of the probes 
were set manually and read against a vernier scale to an 
uncertainty of 0.5°. 

The calibration was carried out for pitch angles 
-15° - $ i 25° in steps of 5° and yaw angles -80° - a - 80°. 150 

measurements were taken while rotating the probe from a = 3 0° to 



41 



a = -80°. During this procedure the probes were yawed con- 
tinuously but slowly enough to acquire data with verified 
accuracy . 

a. Results for the Type A Probe . For an arbitra- 
rily chosen flow condition the output of the probe, p^, is 
shown as a function of yaw angle, a, in Fig. 14 . In Fig. 15 
results for different Mach numbers are given. In this case 
data are plotted only for values - 0 in. H 2 O. In [3] it 
was reported that for a constant pitch angle, (j) , the probe 
reads static pressure (indicated by p^ = 0) at approximately 
the same yaw angle position independent of the Mach number . 
This can be seen in Fig. 15 where the static pressure is 
registered at a yaw angle of a - ±64° for each Mach number. 
This is shown in more detail for an arbitrarily chosen flow 
condition in Fig. 16 , in which a larger ordinate scale has 
been used. At positive yaw angle, all curves cross the 
abscissa at the same location. This is almost but not 
exactly the situation for the negative yaw angle. 

An exception to the preceding statement occurred 
for the lowest Mach number, M = .2, and for the complete 
range of pitch angles investigated. Reynolds number effects 
were suspected, but unfortunately this parameter could not 
be changed independently in the calibration set-up over a 
sufficiently wide range. 

The data reduction procedure for the calibration 
followed, initially, the method described in [3]. For 
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convenience, a FORTRAN computer program was written to carry 
it out. The quantities which must be evaluated are described 
in Section IV. Here is described only how the coefficients 
are determined. 

First, the non-dimensional pressure coefficient, 

c , defined in equation (1), is computed. (Data which are 

P 

obviously incorrect are ignored.) The probe yaw angles, , 
where static (atmospheric) pressure is indicated are evaluated 
by applying a least squares linear curve fit through nine mea- 
sured values which are close to the expected value of . In 
addition to the smallest positive c value the adjacent 4 
positive and 4 negative values are taken. 

A small difference in magnitude was found in the 

positive and negative values of (up to 2 degrees) which 

indicates the presence of a small error in setting the probe 

to zero in the holder when aligned with the flow. Therefore, 

the average magnitude from positive and negative values of 

a is taken, 
s 

The computation of the maximum c^ value, CpQ^/ 
is done by using a least squares curve fit of fourth order 
through the measurement points which are within the range 
-30° ^ a i +30°. When the same method is applied for the 
complete range of a it leads to deviations from the data 
which are largest in the area of the maximum. Simply asking 
for the highest actual value of c can result in a value 
which is too large because of the experimental scatter. 
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(Only obviously wrong values are ignored.) Of course the 



maximum should be expected to occur at = 0. But, again 
because of scatter this is not assumed. Although deviations 
of the value of a corresponding to the calculated value of 
Cpoj^ and the value of a which results from selecting the mid- 
point between positive and negative values of can be up to 
10 degrees, only negligible differences in the value CpQ it- 
self occur because of the flatness of the probe characteristic 
near its area of maximum output. 

As an additional example in Fig. 17 plots for 
Cp/cpQ^ = f(a) are shown for M = . 6 and 25*^ ^ a 1 -10*^. After 
non-dimensionalizing a with , 54 different curves for 
Cp/cpQ^ = f(a/a^) are obtained and a fourth order polynomial 
curve is fitted by the least squares method to each of the 54 
data sets . 



For -1.1 - - +1,1 in steps of .1 2 3 values 

of ^p/^PQj^ calculated for each of the 54 polynomials. 

1242 values are then used to calculate 5 coefficients for a 
single fourth order representative curve. Fig. 13 . Because 
of the fact that the probe behaves somewhat differently at 
M = .2 additionally, coefficients are provided which exclude 
the test at this Mach number. For comparison the result for 
the old A probe is also plotted. 

The deviation of the results for M = .2 is even 
more distinct when CpQ^ is plotted versus pitch angle, <t> , 
Figs. 19 and 20 , and versus ^ , Figs. 21 and 22 . These 
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are additional curves which are used in the reduction of data 
from TX measurements. The data in Figs. 19 and 21 are shown 
connected for each Mach number by straight lines . For each 
of the curves a fourth order curve fit is computed, and with 
the aid of the coefficients 25 values are calculated in the 
range of interest. With the resulting 150 numbers two curves 
as described in the last paragraph are evaluated, Figs. 20 
and 22. 

b. Results for the Type B Probe . For the Type B 
probe only as a function of ip is required [3]. The pro- 

cedure used for acquiring calibration data is the same as was 
described for the A probe. Because of the 35° inclination 
the B probe shows a totally different characteristic in com- 
parison to the A probe. This can be seen from results shown 
in Figs. 23 and 24 , which are representative examples of the 
acquired data. The final result is given in Fig. 25 together 
with the coefficients. The straight lines reqresent the ex- 
perimental data and the dotted ones are the results after 
the curves are fitted. It is to be seen that a distinct Mach 
number dependence occurs. That was not the case for the old 
B probe and, therefore, only one curve is plotted in Fig. 25. 

From a first attempt to calibrate the probe a se- 
verely unsymmetric characteristic was measured which was sub- 
sequently found to be due to some very small dust particles 
in some holes of the screen [45]. When the screen was cleaned 
(which had to be done extremely carefully) a nearly perfect 
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behavior resulted. This experience is mentioned because dust 
particles cannot be avoided completely in any research fa- 
cility like the TX. It is suggested that, when data are sus- 
pect and instrumentation has been checked, this occurrence of 
fouling should be taken into account as a possible cause of 
incorrect measurements. The question can be resolved by a 
visual inspection of the probe tips under a microscope. 

3 . Application of Reduction Programs to Calibration Data 

The programs which were written to evaluate the 
velocity vecotrs from probe measurements were applied first 
to the calibration data. The idea was to provide results 
which would indicate the degree of accuracy in applying the 
approximations. In [3] a similar test is described. There, 
both probes were physically remounted in the calibration rig 
at corresponding yaw and pitch angles and verification of 
different flow conditions v;ere calculated using the reduc- 
tion program. 

In the case of the new probes, remounting and re- 
testing the probes was not necessary since calibration data 
were taken at entirely similar flow conditions for the two 
probes. The reduction procedure was simply applied to the 
data taken during calibrations and the results compared with 
the known conditions set in the test apparatus. 

a. Results for Analytical Data Reduction . Figs . 

26 to 28 show some of the results of several that were ob- 
tained using the method described in Section IV. Diagrams 
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in Figs. 26 and 27 are results from the old probes and those 
in Fig. 28 are from the new pair of probes. Using Eq. (6) , 
for data at fixed Mach number and pitch angle the values for 
Pq^ and Pq 2 are constant and p^(m) and a^(m) are the probe 
pressure and yaw angle. Mach number, M, and pitch angle, 

(p, can be evaluated for each value of Pj^(m) and a^(m) and 
these are shown in Figs. 26 to 28. The fixed test condition 
is indicated as a reference value. It is to be seen that 
the method gives unacceptably large errors in the area of 
aobut -40° < a < +40°, where dotted lines are used to connect 
corresponding results. This is to be expected. But even in 
remaining regions the results are not accurate, especially 
the Mach number which depends (in an iterative process) on 
the evaluation of the static pressure. (The values for total 
pressure are quite acceptable.) In the pitch angle there 
seems to be a better consistency because deviations from the 
test condiitons are reasonably constant in the area where 
good results are expected. 

Those general remarks are true for both pairs of 
probes. For the new probes, much more calibration data are 
available and, therefore, the evaluated coefficients and 
corresponding results were initially expected to be better 
than those for the old probes. But of this there is no in- 
dication in the plots. The reason which is suggested at the 
moment is that the procedure is sensitive to the accuracy of 
the fourth order curve fit. The calibration (vs yaw) is 
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represented quite well for the B probe and for the A probe 
away from the neighborhood of a = 0, but it may be that it is 
not good enough for the A probe near to a = 0 . In [8] a sine 
function was used to represent the calibration data for an A 
probe. The function provided a good fit to the "flattened" 
type of characteristic. It is possible that the use of a 
similar function for the (smaller) new probes would improve 
the results obtained with the analytical method of reduction. 

In Fig. 28 results are shown for reductions car- 
ried out using coefficients evaluated for different Mach num- 
ber ranges. Results for symbols x and + are evaluated using 
coefficients given in Figs. 18, 20, 22 and 25. It is to be 
seen that the curves always show similar trends and are with- 
out significant differences in quality with respect to accur- 
ately reproducing the known test conditions. What is somewhat 
surprising is that even the data reduced using the coefficients 
evaluated for the very same flow Mach number and pitch angle 
condition set in the calibration rig (M = . 5 and (f) = 0°) , 
which are shown plotted with symbol o, show a very similar 
behavior. Therefore, it is likely that the problem lies in 
the particular choice of analytical function and the high 
sensitivity in an analytical approach of this kind. 

b. Results for Numerical Data Reduction . The nu- 
merical reduction method of [14] — Method II, was used to re- 
duce the calibration data for the old A and B probes. The 
available results from the TX were obtained with this set of 
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probes and therefore a test of the procedure for these probes 
was desirable. 

In Figs. 29 and 30 the results for pitch angle, 

, yaw angle, a, and Mach number, M are shown plotted versus 
±a. Using the abscissa in this way should indicate that two 
(of the four required) pressures are taken from the left and 
right branches of the calibration curve for the same absolute 
value of a. The other two pressures are the maximum outputs 
of each of the A and B probes. It should be mentioned that 
every other combination of pressures (which, in principle, it 
is possible to use) gave poorer results in comparison to the 
aforementioned choice. For example, it is possible to use 
three values from the B probe and only one from the A probe, 
and even four values from only one probe will lead to a suc- 
cessful solution of the iteration. It is stated however in 
[14], and now confirmed by experience, that using both maxi- 
mum values of the probes and two others which correspond to 
symmetric yaw angle positions will give the best results. 
Indeed, the yaw position which is recommended in [14] is 
r25°, but for real values of calibration data yaw angles in 
the range of about 50° to 65° must be used. This is quite 
clearly indicated in the diagrams of Figs. 29 and 30. It is 
to be seen that in this range of a the required flow quanti- 
ties show an asymptotic behavior (for simplification the re- 
sults are connected by straight lines) and the results are in 
quite reasonable agreement with the test conditions. In 
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areas near maximum probe outputs (±a < 25°) incomprehensible 
results occur as was the case of the analytical method 
(-40° < a < +40°) . 

The scatter of the results is less in comparison 
to that found using the analytical method (Figs. 26 to 28, 
where the same scales are used). Therefore, for the moment 
the numerical method is recommended as the more successful 
one . 

If ±a values are used which are lower than those 
from the recommended range, yaw angle deviations from the 
test conditions increase very rapidly. The behavior of the 
yaw angle calculated using data at large yaw angles is not 
as obviously asymptotic as are the Mach number and pitch 
angle. This can be interpreted as an indication of the sen- 
sitivity of the reduction procedure. In view of the magni- 
tude of the possible error in a it might be better to derive 
the yaw angle as described for the analytical method in an 
earlier section from the position of the maximum of the probe 
output. It is suggested that these different results for 
yaw angles be compared. This topic will be discussed in 
connection with the results from the TX. 

B. RESULTS FROM THE TX 

Results from high response pressure transducers are ob- 
tained with the A and B probe system in the measurement 
plane between rotor and stator, and with twelve different 
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flush-mounted transducers at locations in the outer casing 
wall. For identification the wall static transducers are 
labeled K6 to K13 a number (see Fig. 5) . The number indi- 
cates the axial location, where 9 is at the rotor entrance 
plane and smaller numbers indicate planes upstream from 
there . 

An example of the data obtained from the Type A and B 
probes using the acquisition procedure described in Section 
III are given in Fig. 31 . Some general remarks can be made. 
Outputs for both probes are shown plotted for one blade pair 
for two different numbers of repetitions (compare Fig. 12). 
[Note that all results for the A and B probes presented in 
this report are taken at the mid-span position (50% of blade 
height) . ] The number of counts in the circumferential direc- 
tion is indicated by n; one passage is divided into 123 
counts. The noise reduction in fluctuations which accom- 
panies increasing the number of repetitions is clearly seen. 
Of course, this would improve if even higher numbers of 
repetitions were used. With the software improvements de- 
scribed in [10], and with the number of repetitions equal 
to ten, the data acquisition time has been reduced to about 
10 secs from the 45 seconds which were needed with the old 
software in [6]. Additionally, it should be mentioned that 
the averaged value indicated in the diagrams by a broken 
line is nearly independent of the number of repetitions 
used. The wakes of the blades are clearly indicated and it 
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is to be seen that they occur almost exactly 128 counts 
apart, which is another indication of the precision of the 
high speed data acquisition. 

In Fig. 32 pictures taken from an oscilloscope are pre- 
sented to demonstrate transducer outputs for the A and B 
probe and for the wall statics, in this particular case K9.5 
and KIO. It is clear from the l/30th sec. exposure that 
much more noise occurs for the probes inserted into the flow 
than for the measurements made at the wall even if the dif- 
ferent scales are taken into account. However, results in 
Fig. 31 also indicate that increasing numbers of repetitions 
will result in well-defined average signals. This would not 
at first be expected from watching the probe outputs on the 
osci lloscope . 

Results for A and 3 probes can now be discussed in more 
detail. 

1. Results for A and B Probes 

In Figs. 33 and 34 , results for different numbers of 
repetitions are shown (in this case for the B probe only) . 

In these diagrams, following the DPDS procedure described in 
Section IV, signals are shown from the B probe when set at 
specific yaw angles, a (for definition see Fig. 3) . It is 
noted that the averaged values for corresponding curves with 
different numbers of repetitions would differ slightly, be- 
cause data were taken on different days. It is to be seen 
that with increasing values for a the basic behavior of the 
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signals within the wakes changes. For example, in the case 
of a = -43° within the wakes extrema are indicated as minima 
and in the case of a = +103°, maxima occur. The difference 
between peak to peak output values decreases for a-values 
which are between these extremes. 

With the aid of Fig. 35 a qualitative explanation 
of this behavior can be given. Here velocity vectors, probe 
positions and curves of "probe output, p, versus angle, a," 
are sketched. (It is of no importance for the following 
qualitative description but it should be mentioned to avoid 
misunderstanding that these curves are actually different 
for different positions of the probe across the blade pas- 
sage.) If the rotor blade wake is a simple wake, the rela- 
tive velocity in the wake, ^ 2 ^, is smaller than that outside 
the wake, W 2 , but the relative flow angle is almost constant 
This can be seen to increase the absolute flow angles seen 
by the probe (u 2 = outside the wake; U 2 q = inside the wake) . 
This results in probe outputs p ^2 ^.nd P^2 d ^°^ ^ ~ 
sketched on the right hand side of Fig. 35, Here, p ^2 ^ ^c2D 
but both have negative values (that means below atmospheric 
pressure), compare the result for a = -43° in Fig. 34. If 
the probe angle, a, is increased both values p ^2 ^c2D 

are "traveling" on curves p = f(a) in the positive a direc- 
tion. For example, in the case of a = +52° both values are 
quite the same, and therefore distinct peak-to-peak differ- 
ences do not occur. For a = 103° ?^2 d ^ ^c2 explains 

the "inverted" signal behavior. 
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It is recommended in [3] that measurements for seven 
different yaw angle positions should be obtained to enable 
an adequate reduction of data from the DPDS technique. The 
first data were obtained (examples are shown in Figs. 33 and 
34) following this suggestion. However, it was subsequently 
found that more angle settings were required to improve the 
data reduction procedure. Both acquisition program [16] 
and reduction program allow measurements at eleven different 
probe settings, and this feature was used in subsequent 
tests . 

Two sets of data from the TX using the old A and B 
probes are shown in Figs. 36 to 39 . They can serve to il- 
lustrate some general statements with respect to the DPDS 
technique. (Following the failure of the A probe as a re- 
sult of handling, further measurements were not possible 
with this system. ) 

Each section of Figs. 36 to 39 shows results for 
both A and B probes for different blade-to-blade positions 
and each represents a discrete value of the probe angle 
setting. For the first data set in Fig. 36 and Fig. 37 the 
number of repetitions was 10 and for the second set in Fig. 
33 and Fig. 39, 40 repetitions were used. 

Comparing the corresponding signals from the two 
data sets, it can be seen that for the first set an apparent 
"phase" shift occurs between the A and B probe results which 
is not the case for the second data set. In fact a change 



54 



was made in the data acquisition software following the first 
data set. It was found that one probe was placed at a dif- 
ferent circumferential location (135° instead of 90°, Fig. 

5) after a bearing failure of the drive turbine was repaired 
and the instrumentation was replaced. Therefore, the number 
of counts separating the two probes which was programmed in 
the acquisition software (to provide measurements with the 
two probes from the same blade pair) did not match the hard- 
ware any more. After correcting the number of counts in the 
program, the signals became aligned correctly, as seen in 
Figs. 38 and 39. 

There is, indeed, no obvious check which would en- 
sure that measurements are obtained from the same blade 
pairs, but the test with the PFG described in Section V.A.l 
indicated clearly that the trigger system worked correctly, 
when correctly programmed. 

The phase-shift of the first set of data could be 
corrected partially by providing an offset of a fixed number 
of counts in the reduction program. 28 counts were needed, 
but the disadvantage remains that signals were taken behind 
different blade pairs. Furthermore, a complete reduction 
over two passages was not possible, unless the results which 
were offset to negative counts were patched to the data at 
the other end to complete the two passages (a procedure 
which may be doubtful) . 
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The corrected software in [10] needed test times of 
about 25 minutes for one complete set of data from 2 probes 
at a fixed radial location, taken at eleven probe settings 
with 40 repetitions. A third of this time was required to 
manually adjust the probe yaw angles. Although the ratio 
of "test time to quality of signal" is acceptable, it may be 
worthwhile to provide remotely actuated probes in the future, 
firstly to decrease the test time and secondly for the con- 
venience and safety of the operator. 

Applying the data reduction (with on-line calibra- 
tion) to the signals obtained from the A and B probes re- 
sults in 2 X 256 different curves of p = f(a) as shown for 
the first data set in Figs. 40 to 45 . The curves shown are 
the result of plotting probe pressure, p, for different cir- 
cumferential locations, n (counts) for the 11 yaw angles, a. 
Measured data points are connected with straight lines for 
clarity. In Figs. 40 to 45 the curves are drawn to a scale 
which is shown once at the bottom for the lowest curve and 
the other curves each offset in sequence by about 3.4 inches 
of water. In Figs. 42 and 45 the results in the area of the 
wake of the first blade is shown in more detail, with each 
curve offset by almost seven inches of water. 

A second set of data were obtained at a highly throt- 
tled flow condition where a higher degree of flow separation 
can be expected to occur in the rotor blading. The same yaw 
angle settings were used. Plots of p = f(a) for these data 
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indicate that especially within the wakes (particularly for 
the A probe) yaw angles greater than 103° were needed in 
order to fuuly describe the right side of the characteristic 
curve. An example is shown in Fig. 46 . This is a diffi- 
culty because the reduction procedure requires adequate data 
on both left and right branches. 

In wake flow areas the probe signal seems to be 
relatively more unsteady when the probes are set to flow 
angles which are close to the mean flow direction. This can 
be seen in the A probe signals shown in Figs. 46a and b . 

For a probe angle which is close to the aligned flow direc- 
tion (a = 52°, Fig. 46a) it is seen that within the wake 
the noise is relatively more pronounced than it is outside 
the wake. For the same numbers of repetitions, at large 
values of probe yaw angles (a = 103°, Fig. 46b) , the noise 
inside the wake flow appears to be very similar to that out- 
side the wake flow and quite similar to that outside the 
wake flow at the probe angle of 52°. 

In view of this there are two possible ways to in- 
crease the quality of the signals without extending the data 
acquisition time. First, the number of repetitions can be 
increased for probe angles which are close to flow-aligned 
conditions and can be decreased for the rest of the angles; 
or second, for all probe angle settings the number of repe- 
titions can be increased when data are acquired in wake flow 
areas and can be decreased outside those areas. The 
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identification of wake flow areas can be made from a display 
of the transducer output on the screen of the oscilloscope. 

At first glance, a comparison of the curves shown 
in Figs. 40 to 46 with those from the calibration, Fig 17 
for example, will lead to the conclusion that the test data 
result in pressure characteristic curves similar to those 
obtained in calibration tests. And it may be considered sur- 
prising that signals like those shown for example in Figs. 

36 and 37 result in quite well-shaped curves, especially 
outside the wakes; moreso, in view of the unsteadiness evi- 
dent in the oscilloscope traces (Fig. 32). Nevertheless, 
going further into the details it turns out that the initial 
statement does not hold true everywhere, especially on the 
ends of the left and right branches and particularly within 
the wake flow areas. It is known from applying the reduction 
methods to the calibration data that pressures must be used 
which are obtained at probe yaw angles which are approxi- 
mately ±60° away from the angle for maximum probe output. 
Techniques are discussed below which may help to overcome 
this very important problem. 

While being aware of the aforementioned problem in 
practice a data reduction procedure is described next which 
is based on the assumption of having a perfect pressure 
characteristic curve. This is done to gain a better under- 
standing and illustrate the initial problems of the DPDS 
technique . 
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From the first data set, the results at n = 31 (or 
n = 3, if the added offset of 28 counts is not included) are 
shown plotted in more detail in Fig . 47 . It is noted that 
it was carefully verified that applying the offset did not 
change the shape of the probe characteristic. 

Using a fourth order curve-fit to the data points 
shown plotted with the symbol • (through which a curve has 
been hand-faired) and calculating corresponding points using 
the fourth order polynomial results in pressure values indi- 
cated by the x symbols. It is seen that in this case the 
polynomial curve fit was a good representation of experi- 
mental results. Next, a line is drawn parallel to the ab- 
scissa and moved with respect to the abscissa to find 
iteratively the pressure p^ , where the length of the line 
between the sides of the curve is equal to 128°. This pro- 
cedure is used because calibration data showed that static 
pressure is sensed by the A probe at about 64° (or ±64°) 
independently of Mach number and pitch angle (Figs. 16 and 
22) . In Fig. 47 values are indicated for p and , a„ and 
the corresponding midpoint, a^, which is regarded as the 
flow yaw angle. 

In addition, it is assumed that the A probe indi- 
cates total pressure, p^, in its maximum output value, which 
is very reasonable as long as pitch angles, $ , are roughly 
in the range | j) j < 12°. In the ideal case, the maximum 

pressure, p , which is evaluated from the fitted curve, 

^ ^max 



59 



should be located at a^, but a lack of symmetry in the charac- 
teristics is known to be present from the calibration tests. 
The difference between and the angle for maximum output is 
about 7 degrees in the example shown. Because of the flat 
behavior in the neighborhood of the curve’s maximum, p (a^) is 
used as total pressure, p. , instead of p . The ratio p /p 
then allows the calculation of the local Mach number, M. 

The assumption that the A probe is not indicating 

static pressure, p^, at = 128° but ±1° from this 

value will lead to deviations in Mach number of about ±2% 

from the previously calculated value. The general trends of 

the calculated Mach number distribution for a set of data 

would not change but the average value would change. The 

same situation will occur if the yaw angle is evaluated at 

p and p__„ is taken to be the stagnation pressure. 
m3.x rTla.X 

For both Mach number, M, and yaw angle, a results of 
the aforementioned procedure for all n = 256 counts are shown 
in Fig. 48 for the first set of data. The results follow a 
credible trend, but the precision will be questioned in the 
following discussion. 

In principle, the same procedure can be applied also 
to the second data set (at reduced flow, throttled compressor 
conditions) . Some mathematical tricks are needed to overcome 
the problem of missing data in wake areas. This has been 
done; however the results are not presented because the ac- 
curacy is seriously questioned. In comparison to the first 
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results a smoother Kach number distribution, higher yaw 
angles and larger wake widths are calculated, all of which 
could be interpreted as the result of a more severe flow 
separation due to the highly throttled test conditions. 

The yaw angle distribution can also be evaluated 
using results from the B-probe . The characteristic of this 
probe is slightly unsymmetric as is the case for the A probe 
and therefore results for v;ill again depend on the pres- 
sure level at which the calculation is made. Because of its 
extreme inclination to the flow even at zero pitch angle 
static pressure is not well defined at a particular a-value 
and must be based on the results from the A probe. It was 
found that a difference of about 10 degrees (dependent on 
which level of pressure was chosen at which to evaluate it) 
occurred in the values of a calculated using the A or the B 
probes. This was ^true for both sets of data. If all uncer- 
tainties in the various methods of evaluating a were con- 
sidered, the difference (which was consistently there) of 
10 degrees could not be explained. However, a misalignment 
between one or both probes and the angle scale on the probe- 
holders is a probable explanation since no verification of 
the angle zero had been made after several months of opera- 
tion. Unfortunately a failure occurred in the A probe as 
a zero check was being prepared. Since the old A and B 
probe set would not be used further, reasons for the de- 
viations in yaw angles were not resolved. 
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Results calculated here for the flow quantities are 
not considered to be satisfactory using either the analytical 
or the numerical data reduction methods, particularly for 
the pitch angle, (}). It was first thought that the discrepan- 
cies in the data-sets described earlier (the presence of an 
offset between A and B probe data, and incomplete data for 
the right side of the pressure characteristic) might be re- 
sponsible. Therefore, parts of another data-set obtained by 
F. Neuhoff (and reduced using a third data reduction proce- 
dure [45]) was used to check the general behavior of the 
characteristic pressure curves derived from the probe out- 
puts. This particular set of data was entirely free of the 
disadvantages mentioned earlier: data were taken behind the 

same blade pair, 40 repetitions were used and the yaw angles 
covered the required range, even in wake areas. It was found 
that the curves of p = f(a) for the third data set were quite 
similar with respect to their basic shape to those shown in 
Fig. 40. Consequently, applying either reduction method to 
this set of data will lead to similarly doubtful results. 

In can be stated therefore that the basic problem in 
the data reduction so far is that the shape of both left and 
right branches of the characteristic pressure curves ob- 
tained from test data do not fully agree with the behavior 
obtained during the calibration of the probes. Since only 
the maximum pressure value is taken from the B probe the 
problem lies only with the A probe. 
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Because of varying pressure levels accompanying the 
varying Mach numbers and pitch angles which are to be deter- 
mined, a comparison of the calibration curves with the appli- 
cation curves on the basis of raw pressure is not useful. 

For a more reliable check the non-dimensional pressure coef- 
ficient, Cp^, given by Eq. (1) , is examined as a function of 
a. For the calibration data, the plots of ^p^^/^pQ^ = f (ot) 
are shown in Fig. 17. Although the curves are quite inde- 
pendent of the different flow conditions, and the value of 
^POA only on pitch angle, t), and is very close to 

unity for 0° < < 15°, it is not possible to generate a 

direct comparison with the test data unless some assumption 
is made concerning the value of p^ or of Mach number. 

The total pressure, p^, can be evaluated reliably 
using the procedure described in Fig. 47 but the static 
pressure, p^ , in Eq. (1) is not known. Therefore the ex- 
pected .Mach number range was used as a parameter and cor- 
responding Pg values were used in the computation of Cp^. 

This procedure was applied to the data from about 70 counts 
of the third data set [45], which were taken at 9 different 
yaw angles, a. 

In Fig. 49 one example is given for data obtained 
outside the wake flow. This flow condition could corres- 
pond, for example, to n = 1 in Fig. 40. The calibration 
curve Cp^ = f{a), plotted on transparent paper, can be used 
to search "by hand" for the particular curve to which it 



53 



best fits. For the example shown this is the curve for 
M = .42. The flow angle, a, can then be "corrected" in areas 
where the two curves do not fit. The resulting values are 
given in the diagram. It can be seen in Fig. 49 that the 
two smallest and the largest values of have to be "cor- 
rected" and that the "correction" for the smallest value is 
the largest. The curve containing the new values of is 
shown, and it is seen that it agrees very well with the cali- 
bration curve. 

In Fig. 50 , as a second example, results are shown 
for data obtained within the wake flow area (for example at 
n = 43 in Fig. 40) . It is clear that there is agreement be- 
tween the test data and the calibration curves at a higher 
Mach number, M = .48, and to "correct" on the left branch 
requires corrections which are about twice as large as were 
required in the first example. For the right branch approxi- 
mately the same magnitude of correction is required. 

This behavior can be explained with reference, again, 
to Fig. 35. The value of the absolute velocity in the wake 
may increase due to the fact that the relative velocity, ^ 2 ^, 
is reduced. (Note that if W 2 j^ 0 , C 2 q"^u.) Moreover, the 
sketches indicate that for the left branch (for example, 
a = -43°) the flow direction to the sensor of the probe for 
°2D ^ more extreme yaw angle than for c^ which 

may explain the necessity for the large corrections. For 
the right branch the flow yaw angle is more extreme for c^ 
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than for ^ 2 ^, but possibly not as extreme as for C 2 ]-) 

a = -43°, which can explain the use of about the same 
correction. 

It should be mentioned again that the general behav- 
ior demonstrated for these examples is typical of the passage 
positions from which they were taken. A similar behavior 
was found for both data sets shown in Figs. 36 to 39. The 
implication is that the data reduction could be carried out 
with data from fewer repetitions (10 and 40, respectively, 
have been used) if an explanation (and a proper correction) 
for the yaw angle problem can be found. Also an "incorrect" 
behavior of probe pressure with yaw angle was observed for 
the B probe measurements. However, a "correction" procedure 
for the B probe requires an additional step, since Cp^ de- 
pends on (}), as can be seen from Fig. 24. 

There are many possible reasons for the observed be- 
havior. What can be said at the moment however is specula- 
tive because detailed investigations have not been carried 
out. For example, the exact position of the probe tip in 
the flow during rotation must be questioned. It is diffi- 
cult to ensure precise mechanical rotation about the center 
of the sensor tip when the probe is yawed in its eccentric 
mount. If the probe tip is slightly long, there is an error 
in position which increases with increasing probe yaw angle 
(in both positive and negative directions) . Different probe 
mounts are used for calibration compared to compressor 
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applications and while the probes are rigid in construction, 
small differences in bending could occur in the two situa- 
tions. It is recalled that there was no opportunity to re- 
calibrate the first set of probes. Although the dimensions 
of the old probes were as small as possible and yet accom- 
modate the transducer, they were possibly too large for their 
intended application. In Fig. 48 the relative size (tip 
diameter) of the old probes is indicated at an arbitrarily 
chosen circumferential location. The further application of 
the smaller probes will give some information concerning 
this problem. 

Errors in tip position on rotation would result in 
an apparent phase shift of the signals obtained at different 
yaw angles. A wake decay due to probe movement to different 
axial distances is demonstrated, for example in [31]. An 
apparent decay is to be seen, for example, in Fig. 38 where 
the positions for minimum output are slightly displaced for 
the A and B probe results for yaw angles of a - -15°, More- 
over, a comparison of the position of minimum output for 
a = -42° and a = -15° for the probe shows a shift of about 
10 counts. However, it must be remembered that in the wake 
all parameters of the flow are changing and the apparent 
shift in minima may be qualitatively correct. 

It is not possible at the moment to give a full ex- 
planation for the behavior of the probes in the compressor 
flow field. If the potential explanations speculated above 
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are eliminated, then the general question of how high re- 
sponse pressure probes behave for different yaw positions in 
unsteady wake flows must be examined. So far, other investi- 
gations on this topic are not known to the author. 

It should be mentioned that this question is also 
particularly relevant when multiple sensor probes are used. 

In such probes, the transducers are placed in various sur- 
faces away from the probe's face. The different transducers 
are differently oriented with respect to the flow direction 
at all times, so that any difference in response which de- 
pends on orientation will be difficult to account for. From 
Figs. 49 and 50 it is evident here that at about |a| > 40*^ 
departure from the calibration behavior begins. Since mul- 
tiple sensor probes are larger in size in comparison to the 
probes described here it is suggested that their unsteady 
response will depart from their steady behavior. 

The numerical data reduction has been applied to the 
results modified as in Fig. 49, since it was shown in Sec- 
tion V that this method was more reliable in comparison to 
the analytical one (Figs. 26 to 30) . The results are indi- 
cated with a X symbol in Fig. 51 . Comparing the curves with 
those evaluated for calibration data (for example Fig. 29), 
it can be seen that approximately the same tendency occurs 
for the flow quantities <)), a and M. The results given for 
yaw angle, a, are here to be interpreted as deviations from 
the angle, , evaluated with the procedure demonstrated in 
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Fig. 47. As before, the best results are obtained in the 
range of about 50*^ ^ I | ^ 65°. At this point, the data re- 
duction has been carried out as usual except that the values 
of a at the right and left branches were corrected. Thus 
the values of c can be calculated from p. and p for every 
data point and those results can be compared with the cali- 
bration curve Cp^ = f(a). When this is done, it is found, 
as might be expected, that deviations occur for about 
|a| < 50°. Since the calibration table covers only the 
range of -70° ^ a - +70° deviations from asymptotic values 
in areas |a| > 65° might also be expected. 

When the method described in [45] is applied to these 
data nearly precisely the same results are found. 

The value obtained for pitch angle, (p , was about 
-2.5°. In [17] it is reported that the time averaged value 
calculated from measurements with the combination probe was 
about (t> = +4.5° and, therefore, results for p within the 
wake must increase considerably to establish an average 
value of about the measured magnitude. The reduction of the 
example given in Fig. 50 showed increases in the pitch angle 
in the wake, but of insufficient magnitude to produce an 
average of 0 = +4.5°. 

This, in turn, raises the question of how reliable 
are the pneumatic measurements, in particular for measuring 
pitch angles. This topic was discussed to some degree in 
Section IV. A distinct problem in the DPDS technique is the 
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way the OLC must be applied to the B probe. In order to find 
out what would happen if the reference pressures of the B 
probe, Pj-g£ g (or equivalently the intercept) were wrong, the 
results were re-evaluated for differences of +2 and -2 inches 
H 2 O respectively. From Fig. 51 it can be seen that this is a 
relatively minor influence on both M and a but it causes a 
change in b of about ±1.5°. In addition, data were reduced 
for Pp^ values which were increased by +2 inches H 2 O when the 
original value of Pg was used. Results were equivalent to 
the case of -2 inches H^O on the B probe and the original 
value of p^^ (+ symbols are used for both cases) . 

The minor influence of these changes on a and M can 
be explained with reference to Fig. 47 although a different 
reduction process is used. A different reference pressure 
on the probes will change only the level of the curve and 
the evaluated yaw angle, a^, would not change. With respect 
to Mach number, an impossibly large change in pressure level 
of say +20 inches H 2 O (p^ = 462.4 in H 2 O and p^ = 427.4 in 
H 2 O) would change M by about 2% only. Such a change of 
total pressure would however result in extremely large 
changes in loss coefficients. 

2 . Results for Wall Static Pressure Transducers 

Very briefly some results are reported which were 
obtained with 12 Kulite pressure transducers mounted flush 
at the outer casing wall. Their peripheral and axial loca- 
tions are indicated in Figs. 1 and 5. Similar measurements 
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from the TX are reported in [5], [12] and [21] where 3 trans- 

ducers were used. The programs to reduce data to contour 
maps used in [12] and [21], respectively, were written in 
BASIC language. A FORTRAN program to include the 12 differ- 
ent signals now available is in preparation. Such a program 
for about 40 transducers is reported in [32]. 

In Fig. 52 a series of signals is shown. The data 
acquisition time was about 10 minutes including measurements 
of all pneumatic pressures and two OLCs . The increase in the 
magnitude of the pressure fluctuations in going from upstream 
towards the rotor in the axial direction is clear, as is the 
rapid decrease at the rotor outlet. Outside the rotor exit 
plane the signal becomes flat and that there are two blades 
is no longer obvious. Signals for K9.5 and KlQ can be com- 
pared with pictures taken from the oscilloscope. Fig. 32, on 
which the time-base can be adjusted to provide a record of 
more than one blade pair. Some differences in signals from 
different passages are evident which are probably due to 
small geometric differences in the hardware. In particular, 
this is certainly the case for the rotor leading edge plane 
at K9 as can be seen from Fig. 53 . The signals from four 
different blade pairs are plotted, and differences in the 
maximoim and minimum values are clearly evident. This is due 
to slight differences in the tip gap and in the chords of 
the blades. It is noted that the tips of the blades pass 
approximately across the center of transducer K9 , on the 
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The largest amplitudes occur for the transducer K9.5, 
Fig. 52. Consequently this location was chosen to locate a 
pneumatic static pressure probe of special shape which was to 
be used in conjunction with the iterative data reduction pro- 
cedure described in [34]. The geometric shape of the probe 
follows exactly the suggestions given in [34]. The probe 
was manufactured so that it could be substituted at any wall 
transducer location. 

A comparison of the measurements obtained with the 
standard pneumatic static pressure tap, the diameter of 
which was 1 mm, with those obtained with the special design 
showed that there was no detectable difference. However, it 
is suggested that for higher rotor speeds with the attendent 
increase in fluctuation amplitudes the comparison should be 
repeated . 

Additionally, signals from one of the transducers 
were acquired with different numbers of repetitions and a 
least squares curve fit method of up to the order of 10 was 
used to approximate each set of results. This was done in 
an attempt to overcome the problem of long acquisition times 
when the initial software [6] was in use. In principle, of 
course, the same procedure could be applxed to A and B 
probe signals. In Fig. 54 results are shown for 5 different 
numbers of repetitions with a curve fit of ninth order. The 
original signal, which is plotted for Rep = 1 only, appears 
as noise about the fitted curve. It is quite remarkable 
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that results for different numbers of repetitions are almost 
the same. The procedure which is a form of data smoothing 
or averaging could be improved either by dividing the results 
for the 256 counts into more than 2 sections or by using a 
more sophisticated mathematical method (which might result in 
capacity problems for the 21-MX computer) . Although the 
software improvements [10] have reduced test time remarkably, 
it may be worthwhile to keep this procedure in mind for fu- 
ture measurements if a speed-up of the data acquisition is 
required . 
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VI. PROPOSAL FOR A FIRST APPLICATION IN THE RGG 



The test facility for Rotating Annular Cascades (RGG) of 
the DFVLR was first used for investigation on cylindrical 
(conical flow angle a = 0°) shaped test-wheels [1 and 2 ], 

In this configuration the test-section, basically, consisted 
of two cylindrical rings forming an annulus. Following this 
program, since 1980 experiments have been underway dealing 
with conical cascade flows. Such flow conditions occur, for 
example, in the last stages of steam turbines, particularly 
in the tip regions. The goal was to obtain information on 
the flow behavior for those rotor areas where it is obvious 
that flow particles can not travel on streamlines which are 
nearly parallel to the machine axis (which is the commonly 
used simplification in flow models) . 

Three new test sections have been built for three dif- 
ferent conical flow conditions: a = 30° (hub and tip), 

a = 45° (hub and tip) and a = 40° (hub) and a = 45° (tip) . 

A sketch of the test section for the second case (a = 45°) 
is given in Fig. 55. It is noted that measurement planes 
and probe settings are chosen in the manner indicated so as 
to be able to make full use of available hardware. Details 
of the complex program which is carried out in cooperation 
with Forschungsvereinigung Verbrennungskraf tmaschinen , e. V. 
(FW) , Frankfurt, West Germany are given (for example) in 
[47] . 
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Basically, the measurement and evaluation method used for 
conical flows is the same one that was developed for cylin- 
drical flows. Of course, somewhat different assumptions are 
involved in the data reduction procedures for each case. In 
the downstream flow field only those quantities which are 
easiest to measure need be evaluated; namely, absolute total 
pressure and absolute total temperature. Assuming that 
there exists a core-flow area in the axial or conical direc- 
tion of the annulus where only minor changes in flow proper- 
ties occur, two dimensional cascade data can be obtained. 
Results which were necessary to prove this assumption for 
the cylindrical flow conditions are given in [1] and [2]. 

In [48] it is demonstrated that for conical flow conditions 
(a = 30°) such a core-flow also exists. Nevertheless, it 
would be highly desirable to have information of the local 
time-dependent flow quantities also. 

Among several motivations the most important is, of 
course, that it would be of more than only academic interest 
to verify that the core-flow areas are fundamentally two di- 
mensional in character, and to use an independent, second 
measurement method to do it. Although in [1] and [2] the 
limits of this core-flow were defined with respect to blade 
shapes, cascade geometries and downstream Mach numbers, it 
is suggested that for different cascades (= test wheels) a 
check of those regions should be done. Of course, it is not 
necessary to measure point by point throughout the different 
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flow fields of interest, but some limited investigations 
should be carried out. It is expected that these areas can 
be evaluated both more accurately and more quickly (with re- 
spect to test time) with high response pressure probes-. 
Therefore, it is suggested that they should be used in paral- 
lel with the existing probes and techniques in the future. 

In addition, conditions can occur where the region of 
nearly two-dimensional flow is quite small. This has hap- 
pened so far for extreme flow conditions (high supersonic 
relative downstream Mach numbers, in particular), in the 
cylindrical case as well as in the conical one. Results 
carried out up to now indicate that the quality of two- 
dimensional flow decreases with increasing conical flow 
angle, a. The investigations for a = 45° are now being 
prepared, and it can be expected that the flow behavior 
will continue this trend. This provides a further reason 
why preparations should be made to make available a second 
measurement method which has the potential of getting more 
detailed information on the flow. 

Finally, it may also be possible to carry out measure- 
ments of general interest concerning the problem of how to 
measure time-averaged pneumatic pressures. The systematic 
investigations which are available so far are limited to 
fluctuation frequencies of about 4 kHz. This is due to the 
relatively small devices which have been used to produce a 
periodically interrupted flow field. The magnitude of the 
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tangential velocities is limiting and the need to have the 
wave length equal to several times the probe dimension has 
placed a practical limit on the available frequency. It 
would be worthwhile to extend this region to ranges of 
about 10 kHz, which are typical characteristic blade pas- 
sing frequencies in many machines. For this purpose wheels 
in the test section of the RGG could serve as nonsteady 
flow generators, only. By using different kinds of probes, 
with a calibrated high response probe as a reference device, 
the comparison of results for the same axial and radial lo- 
cations (of course, in different circumferential positions) 
can be expected to give quantitative results of interest. 
Another advantage of the test-section is that in the down- 
stream flow field, measurements can be made in measurement 
planes which cover a wide range in the axial direction. 

This can not be done in the TX without removing the stator 
and could not be done in many other similar research fa- 
cilities. Doing similar measurements at different axial 
distances could lead to additional information concerning 
the problem of time averaging fluctuating pressures. 

In these investigations it is suggested that an oil- 
filled measurement device such as is described in [33] and 
[34] be included. The reported probe is well understood 
and proven by its originators who have agreed to demonstrate 
its application in the RGG. 
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In the previous section, estimates were made with the 
aid of examples of how an incorrectly evaluated time-averaged 
total pressure can influence results obtained using the DPDS 
technique. Open questions are evident, and so far no other 
experience exists concerning the behavior of the A and B 
probes in areas of high subsonic and supersonic nonsteady 
flows . 

The following remarks concern the application of the 
DPDS technique in the RGG: The first goal should be to 

provide additional qualitative information of the complex 
conical flow field, if possible for the current tests and 
if this can not be done in the time available, for future 
experiments . 

The blade passing frequencies in the RGG are of the same 
order of magnitude as those in the TX. The maximum wheel 
speed is only half (250 revs/sec), but the number of blades 
can be 2.5 times larger. This can occur for a tip-section 
cascade with a pitch-to-chord ratio of 1., low stagger angle 
(about 20°) and low incidence angle (of the same range) . 

This would result in a frequency of f = 12.5 kHz, which can 
be regarded as the highest frequency to be investigated. 
.Assuming f = 10 kHz and 50 blades the maximum speed will be 
200 revs/sec (12000 RPM) which may be a more realistic value. 

For hub section cascades, where lower pitch-to-chord 
ratios are used (about .5) the number of blades is about 100 
and, again, assuming f = 10 kHz will lead to speeds of 100 
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revs/sec (6000 RPM ) , which is about double what has been 
found to be necessary from past experience. 

First of all it is necessary to provide a trigger system 
for the data acquisition. In principal, a variety of solu- 
tions are possible. A sketch in Fig. 56 demonstrates the 
one which is suggested for the RGG. The test-wheels them- 
selves act as timing disks and two photo-electrical sensors 
(NANO-SCAN, Skan-A-Matic Corp., Elbridge, NY) are used to 
provide pulses at one-per-blade and one-per-revolution fre- 
quencies. The installation of the sensors within the inner 
casing (position @ , Fig. 56) is expected to be easy, mainly 
because the sensors have a flexible "snout" which in addi- 
tion can provide an adjustment of the distance between the 
sensors and the wheel. For initial tests, which will be 
carried out without a device similar or equal to the "Pacer," 
only the "one-per-revolution" signal is needed. 

On the back side of the test-wheels (see Fig. 56) radial 
straight lines are required to be coincident with the cen- 
ters of the blades. This must be done accurately to provide 
a stable input to the pacer and permit data acquisition be- 
hind blade pairs of interest independent of minor wheel 
speed variations. Both one-per-blade and one-per-rev. time- 
raark-radii should be chosen to be as large as possible, to 
provide the most distinct electrical pulses. 

The complete arrangement can be checked out using a disk 
which represents the wheel and speeding it up to test 
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conditions (up to 12,000/min) outside the test section. 

Such tests would provide answers concerning the optimum dis- 
tance between sensor and wheel (a first estimate can be taken 
from data given in Skan-A-Matic ’ s handbook) and will also 
determine what width of the marks in circximf erential direc- 
tion is necessary. If the width is small in comparison 
to the distance between two adjacent blades a "wave shaper" 
section in a pacer device (see Section III) may not be neces- 
sary. Consequently, the marks should first of all be made 
by simply laying it out with a metal scribe. It is possible 
that reasonable results could be obtained using a small 
groove filled with some kind of reflecting material. 

It may be possible to use the blades themselves to pro- 
vide the time-marks. In this case the second sensor would 
be installed in position Fig. 56 and from position 

only the "one per revolution" pulse would be needed. This 
procedure should be used if it is more complicated than ex- 
pected to provide well defined time-marks for the case 
as described. Complications for case may result from 

the abrasive material used in the outer casing groove to 
guarantee optimum tip clearance. Perhaps particles of this 
material would disable the sensor. In any event, a "wave- 
shaper" would be necessary because the geometry of the blades 
would produce some kind of a bell-shaped signal. 

The necessary investment is approximately 4000. -DM which 
includes two units and a third one for replacement purposes. 
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In Fig. 57 a sketch of a total pressure probe which can 
be used as an equivalent to the A probe is shown. Although 
it is desirable (or may be necessary), that the flow should 
act directly on the head of the sensor and even small vol- 
umes should be avoided because of resonance problems, geo- 
metric restrictions of the hardware do not allow such a 
design at the moment. In particular, the diameter of the 
holes in the casing where the probes can be inserted into 
the flow field is 6mm and, of course, the maximum value of 
the axial extension of a probe must be slightly less than 
this value. A rearrangement is not possible at the present 
time. Therefore, the proposed modified A probe is a compro- 
mise between what is possible and what is needed, however it 
is expected that it can serve adequately for familiarization 
with the new technique. 

As shown in Fig. 57, the probe will rotate about an axis 
through the probe tip. Because of the design, manufacturing 
problems might occur, particularly since the transducers 
must be handled extremely carefully. In the knee-area (Fig. 
57) a hole (for the transducer) and a slot (for the refer- 
ence tube) are necessary to adjust the sensor in the indi- 
cated position. An adhesive secures the transducer in that 
location and is used for closing the hole and the slot. It 
may be worthwhile to use dummy transducers initially to find 
out whether the installation will work in the suggested way. 
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In Fig. 58 a second probe for the RGG is sketched which 



can serve for measurements of the static pressure (see for 
example Fig. 32) . First of all this device should be used to 
obtain wall-static pressures in the different measurement 
planes where probes can be inserted. In addition, radial 
probe traverses in the annulus are possible using the probe 
because the dimensions are chosen to fit the traverse mechan- 
isms. This may be most worthwhile in the upstream regions 
where the pitch angle, d, is expected to be small but, never- 
theless, it may also be of interest to try it downstream. 

In comparison to a transducer which is flush-mounted in 
the usual way at the casing wall the probe in Fig. 58 takes 
more time to build and the additional expense is justified 
only because of its more varied use. It is, however, easier 
to build than the total pressure probe shown in Fig. 57 and 
does have the advantage that the flow acts directly on the 
sensor . 

In addition, a pneumatic equivalent to the static pres- 
sure probe of Fig. 58 should be manufactured. This will 
allow an on-line calibration as discussed in Section IV, 
and comparisons can be made with a static pressure probe 
which is already available. In the annulus downstream, a 
five— hole cone probe (giving time averaged total pressure, 
yaw angle and pitch angle) would provide a better under- 
standing of the time averaged flow field. 
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The RGG can also be used for the calibration of probes. 

In [1], [2], respectively the calibration of a temperature 

probe was carried out. In this case, the test-wheel was 
replaced by calibration rings which produced a contour over 
the inner cylindrical wall. The calibration rings were cal- 
culated as two-dimensional half-nozzles. Each ring was 
instrumented along the contoured surface with tappings to 
obtain surface pressure measurements to compare with the 
design distributions. Calibrations in the subsonic regime 
were carried out using a single ring producing a simple 
contraction. Different Mach numbers were produced by ad- 
justing the flow rate generated by the radial supply com- 
pressor. For two different supersonic Mach mambers, 1.2 
and 1.4, two different rings were constructed to designs 
which were produced by the method of characteristics. 

In the case of conical test-sections, calibration rings 
are not necessary. The test-wheel is to be replaced by a 
disk which serves to complete the inner contour. With dif- 
ferent compressor speeds and different bypass ratios the 
Mach number can be varied because the annulus is acting as 
a nozzle. In the measurement planes and (^ , Fig. 55, 

even supersonic speeds can be set as can be seen from static 
pressure measurements along the inner and outer contours 
[47]. 

The Reynolds number can be controlled independently of 
the Macy number. This is effected by controlling independently 
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the total pressures and total temperatures in the closed 
cycle. The total pressure is regulated by either sucking 
with a rotary compressor of sliding vane type or by supplying 
air through the bypass-duct. The temperature is controlled 
by the cooler following the radial supply compressor [1 and 
2 ] . 

With the devices and test procedures described in the 
previous paragraphs it should be possible to obtain informa- 
tion on the wake-flow behavior and boundary layer effects in 
the inner and outer casing wall areas. Information could 
be obtained initially simply by recording signals photo- 
graphically from the screen of an oscilloscope, as demon- 
strated in Fig. 32. 

So far, no provisions are discussed for the evaluation 
of pitch angles, p. In [48] calculations based on a 3- 
dimensional time marching method are reported which indicate 
a remarkable streamline shift even when the blade height is 
relatively small. In the example given the flow was de- 
flected from = 45° at the inlet plane to U 2 = 30° down- 
stream. These numbers may be used to estimate the a-ranges 
which may possibly occur in the experimental conditions. 
However, it is recommended that the techniques described 
above be applied first, and a device to enable pitch angle 
measurements should be designed later when a background of 
experience in the RGG has been established. 
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VII. CONCLUSION 



The Dual Probe Digital Sampling (DPDS) Technique for 
measuring three dimensional flow fields in turbomachines has 
been described and its development at the Turbopropulsion 
Laboratory in recent years has been reviewed. Analysis of 
the first results in the transonic compressor at the Naval 
Postgraduate School was reported. 

Two different data reduction procedures--an approximate 
analytical and a numerical one — were examined. When tested 
on data obtained in free jet experiments, the numerical 
method, which used calibration data arrays rather than 
analytical approximations of calibration data, was found to 
give the better results. This was explained as being the 
result of the inability of fourth order polynomial functions 
used in the analytical procedure to adequately represent the 
data over the required range of yaw angles. Nevertheless, 
the application of the numerical method of reduction to data 
from the compressor did not allow the velocity vector dis- 
tribution to be determined with acceptable certainty. 

In the analysis of the initial results, which is reported 
here, the behavior of high response probes as a function of 
yaw angle in a steady flow free jet and in the unsteady flow 
behind the rotor were shown to be different. The differences 
were greatest for data acquired within the clearly unsteady 
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region of the rotor blade wakes. Differences were least for 
data acquired within the nearly periodic region outside the 
blade wakes. Consequently it was concluded that measurements 
obtained in the compressor can not be reduced at the present 
time by a straightforward application of the reported methods 
to the acquired data using the calibration carried out in 
steady flow; to do so will lead to incorrect results with an 
unknown uncertainty. This result is thought to have impor- 
tant implications concerning the interpretation of unsteady 
data from more complex multiple-sensor probe geometries. 

Speculative reasons for the differences in behavior were 
offered; however, time was not available to the author to 
investigate further. The full understanding must be left as 
the goal of the future program. It is suggested that close 
collaboration between and coordination of the programs at 
the Naval Postgraduate School and at the DFVLR-AVA Gottingen 
Fluid Dynamics Branch, where similar measurements are to be 
made in rotating cascades, will lead most efficiently to the 
required understanding and to further development of the DPDS 
and other high response probing techniques. 
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Fiaure 1: Transonic Compressor (TX) 
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Figure 3: Arrangement of probes in the TX 
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as viewed from drive-end (with Q° case angle) : 
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Fi;^ure 5 : Configuration of probes and locations of case-wall transducers 
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Figure 6: DPDS Technique 
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Figure 7: Data Ac'-i 



99 




Figure 8 



Schematic of 



Pacer 



CO«lPUT£lt 



TRANSDUCER 





1 PER REV. 
OPTICAL 



1 PER BLADE 
OPTICAL 



1 PER BLADE-PAIR 
PACER GENERATED 



256 PER BLADE - 
PAIR PACER 
GENERATED 





1 PER REV. 

PACER GENERATED 




Figure 9: Synchronized sampling signals 
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p ^ = 80 in H-,0; f = 4.5 kHz; rep = 10 
‘^ref 2 
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Figure 11; 



Results from the PFG; 
= 80 in H^O; f = 



p , = 146 in H^O; 

^plenum 2 

9.8 kHz; rep = 10 
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Figure 12: Frequency, generated with wave-generator; 



10 knc 
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Figure 13: 



Calibration facility geometry (not ro scale) 
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Figure 14 ; = f(a) 
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Figure 17 ; Cp^/cp^^ versus a; M = .6 
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Figure 20 : CpQ^ = after curve-fit 
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Figure 21 : = f (<!)); experimental results 
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Figure 22 ; 



f((j)); after curve-fit 
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Figure 23 : Cg 

115 



f (a) 




i 



i 



i 



i 



Figure 24 : 



'Pb 



f (a) ; M = .6 
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Figure 25: “ ^(<!>) 
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Figure 26: Results of "analytical data reduction applied to 

calibration data (old A- and 3-probes) 
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Figure 27 : Results of analytical data reduction applied to 

calibration data (old A- and 3-probes) 
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Figure 28: Results of analytical data reduction applied to 

calibration data (new A- and 3-probes) 



120 






121 



21 



r] 




test condition 




Figure 30: Results of numerical data reduction 

applied to calibration data 
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Figure 31 : Results for the A- and B-probe with different repetitions 
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Figure 34; Results for B-probe; Rep = 10 
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Figure 35: Qualitative 
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Figure 36; A and B probe signals; Rep = 10 
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Figure 37: A and B probe signals; Rep = 10 
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Figure 39 : A and B probe signals; Rep = 40 
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Figure 40 : 



= t(a) 



:or A probe (step 1 to 128); Rep = 10 
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Figure 41 : = f(ct) for A probe (step 129 to 256'); Rep - 10 
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Figure 43: = f(a) for B probe ^step 1 to 128); Rep = 10 
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Figure 46a : Results from A probe; a = 52° 
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Ficure 47 ; 



f(a) for A probe (count 31) 
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Figure 48 : M-number and flow angle distribution versus peripheral position 
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Figure 49 : Correction for ct^ values (position: outside wake flow) 
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Figure 50 ; Correction for values (position: within wake flow) 
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Figure 52 ; Case-wall results; Rep = 20 
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Figure 53 ; Signals from K9 
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Figure 54 : Curve-fit for signals from different numbers of repetitions 
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Figure 56 : Proposal of trigger system for RGG 
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Figure 57; Sketch of total pressure probe for RGG (not to scale) 
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